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Qh. Abstract 
D 

*I2 ' We compute the two-loop corrections to the massless quark form factor 7* — > qq and gluon 

form factor H — > gr<7 to all orders in the dimensional regularisation parameter e = (4 — d)/2. 
The two-loop contributions to the form factors are reduced to linear combinations of master 
integrals, which are computed in a closed form, expressed as T-functions and generalised 
hypergeometric functions of unit argument. Using the newly developed HypExp-package, 
these can be expanded to any desired order, yielding Laurent expansions in e. We provide 
expansions of the form factors to order e 2 , as required for ultraviolet renormalisation and 
infrared factorisation of the three-loop form factors. 



1 Introduction 



The infrared pole structure of renormalised multi-loop amplitudes in dimensional regularisation 
with d = 4 — 2e space-time dimensions can be predicted from an infrared factorisation formula, 
which was first conjectured in [1], where it was formulated up to two loops. A proof of the 
formula, together with an explicit formulation up to three loops was derived later in [2]. The 
simplest multi-loop amplitudes where the infrared factorisation formula can be applied are three- 
point functions, involving two partons coupled to an external current: the quark form factor 
7* — * QQ an d the gluon form factor H — > gg. The QCD corrections to these form factors can in 
particular be used to fix a priori unknown constants in the infrared factorisation formula, thus 
enabling an unambiguous prediction for multi-loop amplitudes involving more than two external 
partons. 

In the infrared factorisation formula for a given form factor (or more generally for a given 
multi-leg amplitude) at a certain number of loops, infrared singularity operators act on the 
form factor evaluated with a lower number of loops. The infrared singularity operators contain 
explicit infrared poles 1/e 2 and 1/e. They do therefore project subleading terms in e from the 
lower order form factors. 

At present, two-loop corrections to the massless quark [3] and gluon [4] form factors are 
known to order e°. Two- loop corrections to this order were also obtained for massive quarks [5]. 
The infrared structure of the massless form factors and infrared cancellations with real radiation 
contributions are described in detail in [6]. Very recently, results to order e 2 were obtained for 
the quark form factor [7]. 

The calculation of these corrections proceeds through a reduction [8-11] of all two-loop Feyn- 
man integrals appearing in the form factors to a small set of master integrals. The reduction is 
exact in e, such that the evaluation of the form factors is limited only by the order to which the 
master integrals can be computed. The massless two-loop form factors contain three two-loop 
master integrals, which can be computed either using various analytical methods [12] or numer- 
ically order-by-order in their Laurent expansion using the sector decomposition algorithm [13]. 
Up to now, exact expressions were known only for two of these master integrals, while the third 
(the so-called two-loop crossed triangle graph) was known only as a Laurent expansion up to 
finite terms [14]. 

In this letter, we derive an exact expression for the two-loop crossed triangle graph in terms 
of generalised hyper geometric functions of unit argument in Section |2J Using the HypExp- 
package [15] for the Laurent expansion of generalised hypergeometric functions, this can be 
expanded to any desired order in e. Together with the exact expressions for the one- and two- 
loop quark and gluon form factors in Section |21 this allows the expansion of these form factors to 
higher orders in e. For illustration, we list the one-loop form factors to order e 4 and the two-loop 
form factors to order e 2 in Section^ these orders appear for example in the infrared factorisa- 
tion of the corresponding three-loop from factors. Finally, Section El contains conclusions and 
an outlook. 
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2 Two-loop master integrals 



The virtual two-loop vertex master integrals were first derived to order e° in [14] in the context 
of the calculation of the two- loop quark form factor [3] . All but the crossed triangle graph A 6 
can be expressed in terms of T-functions to all orders in e. 
Factoring out a common 
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No exact expression for A e was known up to now. Following the steps outlined in [14], we obtain 
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While A2 ) LOi ^3 an d ^4 can be expanded using any standard computer algebra programme, 
the expansion of Aq requires the expansion of generalised hyper geometric functions in their 
parameters. For this purpose, a dedicated package, HypExp [15], was developed recently. Using 
this, we obtain the eighth-order expansion: 
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where we encountered a multiple zeta value in the last term. 
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3 Quark and gluon form factors at two loops 

The tree-level quark and gluon form factors are obtained by normalising the corresponding 
tree-level vertex functions to unity: 
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The unrenormalised one-loop and two-loop form factors are calculated from the relevant Feyn- 
man diagrams. Using integration-by-parts [8] and Lorentz invariance [10] identities (which can 
be solved symbolically for massless two- loop vertex integrals, see the appendix of [16]), these 
can be reduced [9-11] to the master integrals listed in Section El 
The unrenormalised one-loop quark and gluon form factors read: 
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where N = 3 is the number of colours and g is the bare QCD coupling parameter. 

The unrenormalised two-loop quark and gluon form factors for Np massless quark flavours 

are: 
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The renormalised form factors are obtained by introducing the renormalised QCD coupling con- 
stant and the renormalised effective coupling of H to the gluon field strength [4] , and subsequent 
expansion in powers of the renormalised coupling. 
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4 Expansion of two-loop form factors 



The renormalised form factors are expanded in the renormalised coupling constant. In the MS 
scheme, the bare coupling ao = 5 ,2 /(4vr) is related to the renormalised coupling a s = a s (/i 2 ), 
evaluated at the renormalisation scale fi 2 by 
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where 

S e = (47r) e e" e7 with the Euler constant 7 = 0.5772 . . . 

and /Xq is the mass parameter introduced in dimensional regularisation to maintain a dimension- 
less coupling in the bare QCD Lagrangian density. For simplicity, we set /j? = q 2 . If the squared 
momentum transfer q 2 is space- like (q 2 < 0), the form factors are real, while they acquire imag- 
inary parts for time- like q 2 . These imaginary parts (and corresponding real parts) arise from 
the e-expansion of 

A( q 2 ) = (-sgn(q 2 )-i0)-t . (13) 



The renormalised form factors can then be written as 
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Expanding the first and second order coefficients of the form factors to e 4 and e 2 respectively, 
we obtain: 
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5 Conclusions and outlook 

In this letter, we computed the two- loop quark and gluon form factors to all orders in the 
dimensional regularisation parameter e. The principal ingredient to this calculation is the two- 
loop crossed triangle graph Aq , for which we computed an exact expression in terms of generalised 
hypergeometric functions of unit argument, which can be expanded to any desired order in e 
using the HypExp-package. 

A potential application of the form factors derived here is the extraction of the complete 
set of infrared pole terms of the genuine three-loop quark form factor from the recently derived 
three-loop splitting and coefficient functions in deep inelastic scattering [17]. In turn, these 
allow to fix the yet unknown hard radiation constants in the infrared factorisation formula 
at three loops. Parts of these constants were derived previously from J\f = 4 supersymmetry 
relations [18, 19]. 

The two-loop vertex master integrals feature as subtopologies in the reduction of the three- 
loop form factor contributions, appearing if one of the three loops is disconnected from the 
others by pinching the connecting propagators. In this case, their terms to e 2 are required. 

The calculation presented here illustrates the applicability of the HypExp-package in the 
calculation of multi-loop corrections in quantum field theory. Functions similar to those which 
were expanded here appear also in multi-particle phase space integrals in massless [20] and 
massive decay processes [21]. Since all these integrals correspond to particular cuts of multi- 
loop two-point functions, one might expect that three-loop and four-loop two-point functions 
could also be expanded using HypExp to high orders in e, as required for multi-loop calculations 
of fully inclusive observables [22]. 

Note added: While finalising this letter, an independent paper addressing very similar is- 
sues appeared. In hep-ph/0507039 [7], Moch, Vermaseren and Vogt compute the two-loop quark 
form factor to order e 2 and apply it in the extraction of the pole parts of the three-loop quark 
form factor from deep inelastic coefficient functions. In this paper, the hard radiation constants 
for infrared factorisation at three-loops and related resummation coefficients are extracted for 
processes involving quarks only. Expanding our unrenormalised quark form factor (|l()jl to order 
e 2 , we confirm the result (B.l) of [7]. 
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